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We consider the thermodynamics of an infinite-range Ising p-spin glass model with an ad- 
ditional r-spin ferromagnetic interaction. For r = 2 there is a continuous transition to a 
ferromagnetic phase, while for r > 2 the transition is first order. We find both glassy and 
non-glassy ferromagnetic phases, with replica symmetry breaking of both the one step and full 
varieties. We obtain new results for the case where r = p > 2, demonstrating the existence of 
a non-glassy ferromagnetic phase, of significance to error-correcting codes. 

1 Introduction 

Competition between quenched randomness, frustration, and a tendency to order uniformly causes 
non-trivial behaviour in many-body systems. Typical examples include spin glasses with ferromag- 
netic bias , neural networks [Q , proteins Q , and error-correcting codes . When the quenched 
randomness dominates, these systems exhibit randomly frozen (spin glass) states. When the or- 
dering tendency dominates, they exhibit uniformly ordered (non-glassy ferromagnetic) states. In 
the intermediate regime they can exhibit mixed (glassy ferromagnetic) states. The phase diagrams 
around this crossover often have a rich structure. 

One of the principal purposes of this paper is to present a new method to construct the phase 
diagrams of such systems from an understanding of the related system with only the random 
term and an external field term. We apply this method to an Ising spin glass with a p-spin 
Gaussian random interaction and an r-spin ordering interaction, both of infinite range. It has 
been shown |0] using replica theory that for p > 2, the system with neither ordering nor field 
terms has two phase transitions: a discontinuous one step replica symmetry breaking, followed by 
a full replica symmetry breaking at lower temperature. Subsequent work investigated the one step 
transition with the addition of either a field |^ or an ordering interaction with r = p . We extend 
the results to the case of general r and also examine the transition from one step to full replica 
symmetry breaking. This work follows on from a previous paper ||l0[| in which we dealt not with 
an Ising spin system, but with a technically simpler spherical spin system. 

The order of the paper is as follows. In §^ we present the general method. In we review the 
replica theory of the Ising p-spin glass. In §^ we discuss the phase diagrams of the Ising p-spin glass 
with r-spin ferromagnetism. In we present some exact results using a gauge transformation for 
the case where r — p. Finally, we make concluding remarks and present technical appendices. 

2 The constrained magnetization approach to mean field 
ferromagnet ism 

We consider a general system of scalar spins with a Hamiltonian 



where the first term is the Hamiltonian of any system (referred to as the base system), and the 
second term is an infinite-range r-spin ferromagnetic interaction. To leading order in 1/N the 
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partition function is 



JdM exp-NP (^/basc(M) - ijoAf 
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(3) 



where 



/base(M) = In ^basc(M) , (4) 

^basc(M) = V 5 ( A/ - ^ V a, ) exp -/3Hbase . (5) 



{cr;} \ i / 



We recognize /basc(Af) as the free energy per site of the base system when the magnetization per 
spin is constrained (via the (5- function in the partition function) to be Af . Given a knowledge of 
this function, we can carry out the integral over M in (||) to leading order in 1/N by the saddle 
point method, giving a free energy per spin of 

friJo) = -^InZriJo) = /basc(Af.) - -JqM^ (6) 

Np r 

where the equilibrium magnetization Mr{Jo) minimizes the right hand side, and so is a solution to 
the equation 

fLseiMr) - JoM^'. (7) 

We note further that if we take r — 1, the second term of (|l|) represents an external magnetic 
field h = Jo, and so the free energy of the constrained base system is related simply to that of the 
unconstrained base system in an external magnetic field by 

/field (/l) = /baso(Affield(/l)) - /lAffield(/l), (8) 
/^,,,(Mfield(M)-/l- (9) 

Therefore a knowledge of the free energy of either system is enough to determine that of the 
ferromagnetic system via (|^) and : the ferromagnetic term acts as an effective field 

h = JoAf-i (10) 

which must be determined self-consistently. 



3 The replica theory of the Ising p-spin glass 

In the remainder of this paper we shall concentrate on the example of an infinite-ranged Ising 
p-spin glass model. The Hamiltonian we use is 

Hr = ^ ^ii...ipCrn • ■ -O-Jp ^' E CTil-'-'^V (11) 

il<i2---<ip ii<i2---<ir 

where Ji-^...i are independent quenched random couplings given by a Gaussian distribution with 
zero mean and variance plJ^ /2NP~^ , and ai = ±1 are Ising spins. For p = 2 we recover a 
generalized Sherrington-Kirkpatrick model ||ll|,|l2| which has only full replica symmetry breaking. 
In the present paper we consider p > 2 so that there also is a region of one step replica symmetry 
breaking. In the limit p oo, we recover the random energy model (REM) which has no 

full replica symmetry breaking: this will be discussed further in §0. 



In the spirit of the previous section, we consider first the simpler 

"Afield = ^ Jii...ip'^ii ■ ■ ■ — h aj . 



(12) 
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Previous studies [@-^ used the rephca method [0,^ to obtain the free energy per site /field given 

by 



/3 /field = hm - 

n^Q n 



a=tb 



(13) 

which must be extremized for the repUca overlaps qab — jj '^ii'^i'^t) foi' 1 ^ a b ^ n. In the 
replica symmetric (RS) ansatz {qab — q for all a ^ b) this is 



^ „(RS) 
A^ifiold 



J e-'-'/^ In 2 cosh (^^J ^pf3W^qP-^ z + /3/i^ (14) 



where the overlap q is given by the self-consistency equation 

dz 



q = 
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/2 tanh^ J\pf3^J^qP-^ z + ph] . 



(15) 



This RS solution is stable against small replica symmetry breaking Almeida-Thouless fluctuations 
II if 

ip(p - l)/32 j g-.V2 secli4 (^^Ipp2j2qp-l^ + < 1 . (16) 

[1 is 



The magnetization of the RS solution is 

dz 



M = 



(17) 



In the one step replica symmetry breaking (IRSB) ansatz ||l6| the free energy per site is /goM^^^ 
given by 



PL 



where 
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f ifL e--?/2 cosh" G (18) 
J V Stt 
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^ ip/3V2gP-i zo + ^ ip/3V2(gr' - qr') + Ph , 



(19) 



and the mutual overlap qo, the self overlap qi, and the break point x (i.e. the weight of qo in the 
averaged overlap distribution function) are given by the self-consistency equations 

2 
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qi = 



dzp 
V2^' 
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f cosh^G tanhG 



f d£k e-^?/2 cosh" G 
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This IRSB solution is stable against further replica symmetry breaking fluctuations if 



1 . r dzn 2, f^e-^?/2cosh"-^G 



Where the IRSB solution becomes unstable, we expect a continuous transition to a solution given 
by the full Parisi ansatz, which we do not consider in detail here. The magnetization of the IRSB 
solution is 

/. „ f^e-^i/^cosh^G tanhG 

M = / ^ . (22) 

J ^ / ^ e--?/2 cosh" G ^ ^ 

4 Full results and phase diagrams 



The phase diagrams for the Hamiltonian (|11|) were found in two stages. They are shown for p = 5 
and r = 1, . . . , 6 in figure ^. 

Firstly, the free energy and magnetization of the system in a field but without ferromagnetism, 
as given by were found. The RS solution was found by numerical solution of ([l5|), the IRSB 
by numerical solution of (^0|). No attempt was made to determine the solution with full replica 
symmetry breaking (FRSB) in the Parisi scheme, but the onset of unstable modes in the IRSB 
solution, as given by (pl|), is taken to signify a continuous transition to FRSB. 



The resulting phase diagram for p = 5 is plotted in figure 1(a). This Hamiltonian is equivalent 



to (11) with p = 5, r = 1, and Jq = h. There are three phases: an RS paramagnet, a IRSB spin 
glass, and an FRSB spin glass. The replica symmetry breaking transition curve has been found 
previously Where the field is less than a critical value, there is a discontinuous one step RSB 
transition (DIRSB) where x = 1, and the transition temperature rises with h. Above the critical 
value, the transition is continuous, first continuous one step (CIRSB) where — qi, and then 
FRSB; and the transition temperature falls with h, vanishing asymptotically as /i — > oo. (Note 
that the DIRSB transition is continuous in the thermodynamic sense.) The numerical calculation 
of the IRSB/FRSB curve is troublesome near the CIRSB transition line: this is because the 
difference between the left and right hand sides of (21) is very small, of the order of {qi — go)'*- 



Consequently we do not have reliable results in this region, and a small gap is left in the plot to 
indicate this. However, the point where this curve meets the CIRSB line (which we shall call the 
RSB triple point, since RS, IRSB, and FRSB phases meet at it) has been determined accurately 



using the RS solution and a perturbative calculation, details of which are given in Appendix A. 
The magnetization M was calculated with (|l^) and (p2|), and curves of /i as a function of M at 
fixed temperature for p — 5 are shown in figure ^. We indicate where the solution is of each type 
(RS, IRSB, FRSB), but note that h is not exact in the FRSB region, since it is approximated by 
the IRSB value. 

Secondly, the methods of §|| were used to calculate the transitions in the full Hamiltonian, 
essentially by finding the stable solutions of ( [Tc| ) with the h(M) found previously. (In the FRSB 
region, the use of the IRSB values introduces a small error in the ferromagnetic transition lines.) 



The resulting phase diagrams for p = 5 and r = 2, . . . , 6 are plotted in figures l(b)-l(f). 

For r = 2 and p > 2, the phase diagram contains six regions. For small Jq, there are three 
M = phases: on cooling, the RS paramagnet undergoes first a DIRSB transition, and then an 
FRSB (not shown explicitly in figure]^). Increasing Jq, the Af = phases become unstable, and 
a continuous transition to M ^ occurs. There are three ferromagnetic phases: RS, IRSB, and 
FRSB. Within the ferromagnetic region, the RS/IRSB curve has DIRSB and CIRSB parts, as for 
■r = 1. (The RS/IRSB transition in the M = Q region is always DIRSB, since the effective field 
JoAf^-'^ vanishes.) 

For r > 2, the ferromagnetic transition is first order: on increasing Jo, there is first a spinodal 
transition, at which a metastable M ^ Q solution appears, followed by a thermodynamic transition, 
at which the comparative free energies make it the favoured state. 

Increasing r at fixed p, we identify two qualitative changes, illustrated in figure 0. For r 



slightly greater than 2, the diagram has the form of figure 3(a) . On increasing r, the peak of 
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Table 1: Special values of r for various choices of p. There is no DIRSB curve in the ferromagnet 
for r ^ ri , and no IRSB region in the ferromagnet for r ^ r2 . 



the RS/IRSB curve approaches the spinodal ferromagnetic transition hne, and the DIRSB curve 
between them shortens. The peak reaches the spinodal line for r = ri, and there is no DIRSB 
curve in the ferromagnetic region for r ^ ri, as in figure [3(b)| . The RSB triple point approaches 
the spinodal ferromagnetic transition line on increasing r, and the IRSB region shrinks. The triple 
point reaches the spinodal line at r = r2, and there is no ferromagnetic IRSB region for r ^ r2, 
as in figure 3(c) . These special values of r may be evaluated for given p with some accuracy, as 
discussed in Appendix B . Some results are shown in table |^. They are in accordance with an exact 
result, to be proved in §|5|, that ri cannot exceed p. 

Of particular interest is the behaviour of the model in the case r — p: this is relevant to error- 
correcting codes and is amenable to an exact calculation on the Nishimori line. We observe 
that there is a first order transition to a ferromagnetism which is glassy for small enough values of 
Jo and T, and that there is a triple point where the spinodal ferromagnetic transition line and the 
glassy transition line within the ferromagnet meet with infinite slope. (Exact results in the r = p 
case, including the fact that the solution is RS on the line p/SJ^ = 2 Jq, are given in §||.) The RSB 
region shrinks as r = p increases. The existence of an RS ferromagnetic phase persists in the limit 
r = p ^ oo. Our results therefore contradict conclusions of Dorlas and Wedagedera ||l^, who claim 
that the ferromagnetic region in the model with r = p —^ oo is entirely glassy. This discrepancy 
is due to the order in which the limits r p and p ^ oo are taken. The approach of Dorlas and 
Wedagedera is to consider the REM with general r, effectively taking the limit p ^ oo in the Ising 
model, and to subsequently take the limit r ~f oo. In doing this, they obtain results which would 
apply to p ^ r — !■ oo. In our approach we set r = p and then take p —^ oo, obtaining results 
which apply to r = p oo, the appropriate case when one considers error-correcting codes 



This difference may be illustrated using figures 1(c) and 1(e): the former is for r < p, and the 
ferromagnetic phase is dominated by the IRSB region in the limit p — > oo with r < p; the latter is 
for r = p, and the ferromagnetic phase is dominated by the RS region in the limit r ^ p oo. 



5 Special features on the Nishimori Une 

Some of our conclusions about the phase diagram for r = p are reinforced from a different point of 
view. We show in the present section that on the Nishimori line, the distribution function of the 
spin glass order parameter consists only of two delta functions. This immediately implies that the 
ferromagnetic phase is not glassy on the Nishimori line for any r = p. 

We start our argument from the definition of the distribution function of magnetization under 
the Hamiltonian ( 11 ) with r — p: 



E{..}%-A^"'E,^0exp(-/37^p) 



(23) 



E{^,}exp(-/37^p) 

Here the outer brackets denote the configurational average over the Gaussian distribution of ex- 
change interactions. To define the distribution function of the spin glass order parameter, it is 
convenient to introduce two replicas of the same system, with spins {crp"*} and {cp"*}: 



E{<,a,j Hy - ) exp (-^Hp({af )}) - /3Hp({af )})) 



exp (-/^^^({af )}) - /37^p({af ^j) 

Our main result is the identity 

Prniy) = Pq{y) whcrcvcr p/3J^ = 2 Jo. 



(24) 
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Figure 1 : The phase diagrams calculated in §^ for the systems given by (a) the Hamiltonian ( p^ ) 
with p = 5, and (b)-(f) the Hamiltonian (|ll|) with p — 5 and r = 2, . . . , 6. The phases are labelled 
by their symmetry breaking (RS, IRSB, and FRSB) and, for r 2, by their ferromagnetism 
{M = for the paramagnetic (PM) and spin glass (SG) phases, while M ^ for the ferromagnetic 
(FM) phases). For r > 2, where the ferromagnetic transition is first order, the spinodal transition 
is shown with a solid line, and the thermodynamic transition with a dashed line; pha ses labelled in 
italics are metastable. The RSB triple point (found by the perturbative calculation of Appendix A ) 
is shown as a circle. (The transition to FRSB is at a temperature too low to appear for small h or 
Jo, in particular within the spin glass phase where the effective field vanishes.) 




Figure 2: The field /i as a function of the magnetization M for the system given by the Hamiltonian 
(|l|) with p = 5, as calculated in at various temperatures. The solid curves show RS solutions, 
the dot-dashed curves IRSB, and the dashed curves FRSB; the points mark the transitions. The 
field here is identical to the effective field JqM^~^ in the system with Hamiltonian dO), as discussed 
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Figure 3: Phase diagrams for the system given by the Hamiltonian (|ll|), illustrating the qualitative 
changes on increasing r at fixed p. The ferromagnetic phases are labelled according to their replica 
symmetry breaking. For clarity, the paramagnetic and spin glass phases are amalgamated and 
the RSB transitions in the M = solution are not shown. The spinodal transition and the RSB 
transitions in the ferromagnetic phase are labelled in italics. The peak in the RS AT line and the 
RSB triple point are also labelled where present. On increasing r, both these points move towards 
the spinodal line. 
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Figure 4: A sectio n of the phase diagram for the system given by the Hamiltonian (|11| ) with 
r = p = 5, as figure 1(e) , showing the Nishimori line p/3J^ = 2 Jq as a dot-dashed line. The solution 
is replica symmetric with M = q on this curve, as discussed in §^ The lower-right point marks 
where the paramagnetic, spin glass, and RS ferromagnetic phases meet in thermodynamics. The 
upper-left point marks where the paramagnetic, the RS ferromagnetic, and the FRSB ferromagnetic 
phases meet in spinodals. The Nishimori line passes through both points. 



The condition p(3J^ = 2Jo defines a curve in the phase diagram called the Nishimori line [|T^. This 
line passes through the paramagnetic and ferromagnetic phases as shown by the dot-dashed curve 
in figure |4| The result is proved in Appendix C.l. 



It is well established that the distribution of magnetization Pm{v) has a simple structure with 
two delta functions at symmetric positions. Then ( p5| ) proves that the distribution function of the 
spin glass order parameter Pq{y) is also simple on the Nishimori line, which is not the case if the 
ferromagnetic phase is glassy. 

In Appendix C.2 we demonstrate three further facts about the Nishimori line: it passes through 
the point where the spin glass transition line meets the thermodynamic ferromagnetic transition 
line (the lower-right point marked in figure^); at this point, the qi of the spin glass and the q of 
the ferromagnet are equal; and it passes through the point where the continuous glassy transition 
line (or, equivalently, the RS Almeida-Thouless line) within the ferromagnetic phase meets the 
spinodal ferromagnetic transition line (the upper- left point marked in figure ^). This last result 
demonstrates that there is no DIRSB transition line within the ferromagnetic region aX r — as 
the continuous glassy transition reaches the spinodal line which marks the edge of the ferromagnetic 
phases. It follows that ri ^ p, as claimed in §0. 



6 Conclusions 

In this paper, we have shown that the behaviour of a spin glass with an r-spin ferromagnetic 
interaction of infinite range is simple to obtain from a knowledge of the free energy either of the 
base spin glass with a constrained magnetization, or equivalently of the base spin glass in an 
external field (although these latter problems are often numerically subtle); and we have applied 
this method to an Ising glass with a p-spin Gaussian random interaction. We have found that, 
for p > 2 and r ^ 2, there are both glassy and non-glassy ferromagnetic phases; and that for r 
less than a p-dependent critical value the glassy ferromagnet has regions of both one step and full 
replica symmetry breaking, while for higher r it has only full replica symmetry breaking. We have 
also proved several exact results for the case where r — p, which is relevant to error-correcting 
codes: in particular, we have explicitly shown the existence of a curve of simple form in the phase 
diagram, passing through both paramagnetic and ferromagnetic phases, on which the solution is 
certainly non-glassy. We hope that the present analysis will stimulate investigation of the phase 



structures of related systems. 



Appendix A Results from perturbation theory close to the 
Almeida— Thouless line 

In this appendix we derive expression for x^, the value of x along the CIRSB curve, and the 
condition to determine the RSB triple point where the RS, IRSB, and FRSB ferromagnetic phases 
meet. We use the definitions 



Tn = 
Sn = 



-^e-^o/2tanh'' 



/2 sech" ( ^\pf3^J^ql-\.a + Ph 



(26) 
(27) 



A continuous one step replica symmetry breaking (CIRSB) transition occurs on the line where 
9o = 91 in the IRSB ansatz. Details of the IRSB solution close to this line, using only the numerical 
solution to the RS equation, are obtained by an expansion of these equations in small e = (qi — go)- 
Since the hyperbolic functions of G become polynomials in zi , the integrals over zi can be carried 
out immediately, simplifying matters considerably. In particular, one can find the RSB triple point 
using only the numerical RS data. 

To leading order, both (20a) and (EOq) are 



90 = ^2 + 0(e) (28) 

and so go = 9 + 0(e) where q is the RS order parameter satisfying (p^. To obtain a second piece 
of information from these two equations, [( |20b| ) — ( ^Oal )] must be expanded to first order, giving 

e=\p{p-l)l3^j\l-^Sie + 0{ef. (29) 

This shows that (|l|) holds as an equality to leading order in e, and hence that the CIRSB transition 
and the onset of RS AT instability coincide. This is to be expected: where the replica symmetry 
breaking mode of the RS solution softens, a continuous transition to an RSB solution occurs. 

qi. Expanding to first order gives (|28) multiplied by fci 



(20c) is trivially satisfied for go 



\p{j)— l)(3'^J^qQ e, which gives no new information. Expanding [(20c) — fci(20a)] to second order 



gives (pq) multiplied by /c2 = |p(p^ ^){P^ 2)/3^ J^gg plus (]29|) multiplied by ifci, which again 



gives no new information. Expanding [(pOq) — (fci + A;2)(20a) — ^fci((20b) — (20a))] to third order 



gives 



12 



where 



X = Xc + 0{e) 



2{p - 2)g3 + 4p{p - l)/3V2gg+2 - 2[p{p - l)/3^ J^qP]^Se 



[pip 



l)/32j2gP]256 



(30) 
(31) 

(32) 



This gives a perturbative expression for x on the CIRSB line using only the RS value of gj^ 

The calculation of the IRSB stability of the CIRSB transition line follows a similar method. 
Expanding (|2^) in powers of e, everything up to third order vanishes by the self-consistent equations. 
Taking a linear combination of the inequality (|2l| ) and the equations (^) such that all terms to 
third order cancel, the fourth order term gives 

4:{p - 2){8p - I5)q + A8p{p ~ l){p - 2)(3^ J^qP 

- Ap\p - 1)2/34 j4g2(f-i) [4g + 3(p - 2)(2 + x,)S^] 

+ 3p3(p - 1)3/36 j6g3(P-i) [8(1 + x,)S& - (6 + %x, + xl)Sfi] > . (33) 



^Our ( p2[ ) for Xc differs from that given by (36) of de Oliveira and Fontanari [H], due to tfieir omission of one of 
these higher order corrections. Further details may be found in our comment on that paper 



At high temperatures, this is satisfied and the continuous rephca symmetry breaking transition is 
to IRSB. Below a point (called the RSB triple point) found by solving (p^), and ([l^) and ( p3| ) as 
equalities, the transition is to FRSB. Within the RSB phase there is a line marking a transition 
between IRSB and FRSB which meets the RS AT line at this triple point. 



Appendix B An equation for the spinodal transition line 



In this appendix we use the definition 

t = tanh ( ,J^pf3^J^qP-^z + fih 



(34) 



The spinodal transition to ferromagnetism occurs at the minimum Jq for which an i\f 7^ 
solution is possible. In the RS region, the ferromagnetic solution is given by (|l5|), (|l7|), and (10). 
Differentiating these three with respect to /i, and using integration by parts, 



dh 
dM 
'dh 



dz 



V27r 
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-t(l - t2)ip(p - 1)/3VV-'^ + (1 - t')/3 



Asserting the spinodal condition dJo/dh — and eliminating the unwanted derivatives gives 



(35) 
(36) 
(37) 



M ={r- l)(3h 



dz 



2tt 



(38) 



We are often interested in the intersection between the spinodal line and the CIRSB line. At 
this point, (|l^) holds as an equality, and (^8|) simplifies to 



M = (r - l)ph 



dz 



-z'/2 42 



^e-^V2 i2(l _i2) 



(39) 



This result can be used as discussed in §^ The peak of the IRSB transition in the ferromagnetic 
region is given by the solution of the RS self-consistency equation (jl^), the RS AT condition ( p^ ) 
as an equality, and the DIRSB equation = 1 with Xc given by (|3^). When this peak hits 
the spinodal line, the DIRSB curve in the ferromagnetic region vanishes. This occurs at r = ri, 
obtained for given p by solving these three equations and (|39|). Similarly, the RSB triple point is 
given by (p^), (16) as an equality, and the perturbative IRSB AT condition (p3[) as an equality. 



When this point hits the spinodal line, the IRSB ferromagnetic region vanishes. This occurs at 
r = r2, obtained by solving these three equations and (39). Some values are shown in table |^. 



Appendix C Proofs for §^ 



Appendix C.l Proof from gauge transformation 

In this subsection, we prove ( ^5| ) using a gauge transformation technique. We exphcitly write the 
configurational average appearing in (p3|): 
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il<...'jp 

Mp-i)\ 



NP- 



(41) 



where we have introduced a factor of unity at the end, defined DJ = Y[ii< i \j ^pj^ dJi-^ ,,,i^, and 
Zj{f3) is the partition function for a given choice of the couphngs Ji^...ij,'- 



-p 2^ I Jii...if, — — I . . . (J,.. 



NP- 



(42) 



Shifting the integration variable by Jo(p — \)\/Np ^ and reordering trivially, we obtain 



^m(y) = E / ^'^P - 



NP- 



2Jn 



ii<...ip ii<...ip ii<...ip 



E E „72'^*i-«P E 



X exp 



E ^ - ^ E '^^^ exp 
/3 J^i-^p'^n---^: 



-13 Y Jn-.-ipC^'H ■ --^i. 



il<...ip 



(43) 



We now perform a gauge transformation with ai GiTi and Ji^...ij, Jii...ipTi-^ • • ■ ^ip, and obtain 



Pn.{y) = ^ /"pJ exp - 



E '-^n-.-ip E „ 72 •^^i - ^p'^^i ■ ■ • E 



pj2 



Ti exp 



-/3 ^ Ji^.A^'Jii ■■■CTi^ 



X exp 



-/3 E 



(44) 



If 2JqIpJ^ — (3 then exchanging two terms gives 



Praiy) = ^ /PJ exp - 



N'P- 



5Z '^ii-.-ip 



2Jo . 



^ Jo^(p~l)! 

2l<...l„ 



X exp 



/3 ^ ^ Jii...ipTi-i ■ ■ ■ Ti^ 



exp 



-/5 ^ Jit...ipCri^ ■■■CTip 



ti<...ip 



ii<...ip 



(45) 



Finally, shifting the integration variable back, we obtain 



I n PJiJn...^p)dJ^,. 



which completes the proof. 



1 



= PM (46) 



Appendix C.2 Other proofs concerning the Nishimori hne 

In this subsection, we prove the results concerning the Nishimori line and two triple points men- 
tioned in § |5|. (The demonstration of several mathematical identities used here is left to Ap- 



[pcndix C^ .) We first notice that (|5|) means not only that the ferromagnet on the Nishimori line 
(p/3J^ = 2 Jo) is rephca symmetric, but that it has q — M. It follows that \pl3'^ J'^qP~^ = f3JoM^~^, 
and with ( p^ this is /Sh. We may verify q = M independently of the gauge result, as using ( [l5| ) 
and (O) this becomes 



dz 



^ e-^'/2 tanh^ (^^/(3hz + (3h^ ^ J e^^'/^ ^^^^^ (^y^z + , 



(47) 



an identity proved in Appendix C.3 



The spin glass transition line satisfies (^0|) with x ~ 1 and qo = h ^ (as M — 0). Matters 
are simplified by the independence of G from zq, which means that the integrals over zq become 
trivial. The assertio n th at the qi here is the same as the q of the RS ferromagnet on the Nishimori 
line becomes, with (20b), 



dz 



e^^'/^tanh^ (^^/pKz + Phj 



J ^ cosh {y/l3hzi) tanh^ {^/Phzl) 

/^e-^?/2 cosh(v^zi) 



(48) 



which is also proved in Appendix C.2 . The other non-trivial self-consistent equation, ( 20c), becomes 
i(p-l)/3W--ln j -^e-^i/2^osh(y^zi) 



/ e-^?/2 cosh (V^zi) In cosh {^/Phzi) 
/-^e-"?/2 cosh(V;g7^zi) 



(49) 



The thermodynamic transition line is located by equating the free energy of the ferromagnetic 
solution and that of the paramagnetic or spin glass solution. The free energy for r = p is related to 
that in a field via (^ and (H). Where the thermodynamic and spin glass transition lines cross, the 
free energy of the spin glass is /sg = /ficS^^^ with x = I, qo = h = Q, and given by (|l8|). 



If the Nishimori line passes through this point then the previous results apply. The free energy of 
the ferromagnet is therefore 



/] 



FM 



(RS) 



(50) 



with f^^^ given by (|14D. The thermodynamic transition condition /fm = /sG becomes 
i(p-l)/3W- / -^e-'"/Hncosh(^z + ph) = - In / e^^'/^ cosh ( 



or, using @, 

/ cosh (^/^Zl) In cosh {^/pKzi) 

f j|l e-^?/2 cosh (V^zi) 



(51) 



e-^'/2 In cosh (^V/^^-^ + /3/i) = , (52) 



which is also proved in Appendix C.3. Thus, the Nishimori line does pass through the thermody- 
namic triple point where the paramagnet, spin glass, and RS ferromagnets meet. 

The equations for the intersection of spinodal line and the CIRSB transition line are discussed 
in Appendix B| . If this point lies on the Nishimori line then the above results apply. Setting p = r, 
M = q, and ^pfi^J^qP-^ = (3h in (||) and in as an equality, dividing the former by the 
latter, and dividing the result by \p{p — l)/3^J^(j^ , gives a condition for the consistency of this 
hypothesis which may be written 



dz 



27r 



sech^ [yi3hz + (3h^-l + j e"^'/^ tanh^ [y/phz + I3h^ 



J e-^'/2 tanh {^/pEz + (3h) sech^ {y/Wz + Ph) 



J ^ e-^V2 tanh^ {y/l3hz + (ih) sech^ [^z + f3h) 



0, (53) 



which is proved in [Appendix C.3| . Thus, the Nishimori line also passes through the spinodal triple 
point where the paramagnet, RS ferromagnet, and RSB ferromagnet meet. 

Appendix C.3 Proofs of some mathematical identities 



In this section we prove several identities used in Appendix C.2. These proofs will all rely on the 
lemmas 



dz 
dz 



e-^'/2 cosh(az)/(z) = e'''/^ 
e-^'/2 sinh(az)/(z) = e'^'/^ 



dz 



/2tt 
dz 



e-'-'/'foiz + a) 



(54) 
(55) 



where /e and /o are the even and odd parts of the function /. The proof of these runs as follows 
dz 



cosh(az) 
sinh(az) 



f{z) = ll (e''^ ± ) /(z) 

^ L 
/2^ V 



-(z-a)V2 



(z+a)V2 



!{z) 



(56) 
(57) 



ie«V2 
2 ' 



.aV2 



dz\ 



e-'"^^fiz,+a)± 



dz2 
V2^ 



e-'-"-/^fi-Z2 - a) 



dz 



,--V2 



/e(z + a) 
fo{z + a) 



(58) 
(59) 



The proof of now runs as follows: 

dz 



dz 



tanli^(az + a^) = eT" / ^= e^'- 1"^ cosli(az) tanli^(az) using (|5|) 



/27r 
dz 



g-a /2 / ^-z II sijji^j-Q^) tanli(az) 



dz 



e ^ tanh(a2 + a^) 



The proof of (^8|) runs as follows: 

dz _^2,2 , , 
e ' cosh(az) 



- X/2 



/27r 



using (psh 



using 



/27r 



e"""'/^ cosh(az) tanh^(az) = e"'/^ / e "^'^^ tanh^(az + c?) ; using (||) 



J ^^e"^^/^ cosh(az) tanh^(az) 



J d^g_.2/2 cosh(az) 



dz 



e ""^/^ tanh^ (az + a^) 



and the proof of (p2|) is identical with tanh replaced by In cosh. 

cosh(az) sech^ (az) , 



The proof of ( p3[) runs as follows: using (|54), 
/ 4Le-V^sech^(«z + a^) =e-^/2 / ^ 

/ 4Le-^V2tanh^(„, + „2)^g-aV2 / ^ 
using dH), 

--2/2 



cosh(az) , 



= -^V2 



cosh(az) tanh^(az) ; 



e ^ ' tanh(az + a ) sech (az + a ) 



dz 
dz 

7^ 



,-zV2 



sinh(az) tanh(az) sech^ (az) 
cosh(az) tanh^ (az) sech^ (az) ; 



and using again 



dz 



e"""^/^ tanh^(az + a^) sech^ (az + a^) 



dz 



cosh(az) tanh^ (az) sech^ (az) ; 



thus 



dz 



2tt 



dz 



sech^(az + a^) - 1 + / ^= e"^ tanh^(az + a^) 



27r 



+ 



/ -j= e~"^'/2 tanh(az + a^) sech^ (az + a^) 



-aV2 



dz 
x/2^ 



f e-^V2 tanh^faz + a^) sech^faz + a^) 
cosh(az) [sech"* (az) — 1 + tanh^ (az) + tanh^(az) sech^ (az)] ; 



and this is identically zero since the term in square brackets vanishes using sech^(az) 
tanh^(oz). 
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